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Anisotropic  Ginzburg-Landau  theory  for  coupled  s-wave 
and  d-wave  order  parameters  is  used  to  analyze  the 
unique  thermodynamic  and  transport  properties  of  the 
new  La_  x(Ba,Sr)xCu04  and  YBa-Cu^O^  superconductors. 
This  simple  phenomenological  approach  is  used  to 
explain  the  prevalence  of  the  large  Sommerfeld 
coefficients  of  the  specific  heat,  the  existence  of 
multiple  specific  heat  anomalies,  the  ultrasonic 
attenuation  peak,  and  mod^l  the  anisotropic  critical 
field  data  as  observed  in  oriented  samples. 


Following  the  discovery  by  Bednorz  and  Mdller  [1]  of  "high- 
temperature"  superconductivity  in  the  rare-earth  copper  oxides, 
there  have  been  numerous  investigations  of  the  anisotropic 
electronic  [2,3]  and  magnetic  [3-7]  properties  of  these  materials. 
It  is  now  well  recognized  that  any  successful  theory  of 
superconductivity  for  the  high-T  oxides  must  include  the  quasi- two- 
dimensional  nature  of  the  Cu-0  pfanes;  the  theory  must  provide  *  in 
addition,  for  a  coupling  between  the  planes  [8,9].  One  of  the  best 
known  theories  of  the  new  superconductors  is  the  resonating -valence - 
bond  (RVB)  model  of  Anderson  [10]  which  describes  the  onset  of 
superconductivity  as  a  Bose  condensation  of  quasi -particle  pairs 
within  a  large-U  Hubbard  model.  It  has  been  shown  by  Kotliar  [11] 
and  Inui,  et  al  [12]  that  the  superconducting  order  parameter  of 
this  model  possesses  s-wave  and  d-wave  components,  the  latter  being 
favored  at  large  U  and  near  half- filling.  At  low  temperatures  the 
mixed  (s+d) -state  is  favored,  similar  to  that  found  in  the  heavy- 
fermion  superconductor  U.  Th  Be^~  [13-15].  It  is  interesting  to 
note  that  the  low- temperature 1  behavior  of  the  penetration  depth, 
A(T)  [16],  the  large  Sommerfeld  coefficients  of  the  specific  heat,  7 
[17,18],  the  enhancement  of  the  sound  velocity  and  ultrasonic 
attenuation  [19,20],  and  the  thermopowers  [17]  of  the 
La  (Sr.Ba)  Cu04  (called  214)  and  YBa^UjOy  g  (called  123) 


materials  are  very  similar  to  the  heavy- fermion  systems.  This  leads 
us  to  believe,  as  has  been  suggested  on  the  basis  of  high-resolution 
X-ray  scattering  experiments,  [21)  that  s-  and  d-wave  coupling  may 
exist  in  the  high-Tc  superconductors. 

Model 


In  this  work  we  apply  anisotropic  Ginzburg- Landau  (GL)  theory  [22), 
previously  extended  by  us  to  include  coupled  s-wave  and  d-wave 
superconducting  order  parameters  [23] ,  to  qualitatively  analyze  the 
single-crystal  and  oriented-f ilm  data  on  the  214-  and  123 -materials . 
In  particular  we  think  that  the  large„Sommerfeld  coefficients  7  -  5 
mJ/mol  KZ  [4,24,25]  and  9  mJ/mol  K  [18,20]  for  the  40  K  and  90  K 
superconductors,  respectively,  the  anomalous  peak  in  the  ultrasonic 
attenuation  at  T  -  0.9  T  [19,20],  the  upturn  in  the  Hc2(T)  curve 
[6,7],  and  the  anisotropy  Cin  the  magnetic  properties  of  these 
materials  can  be  explained  in  the  context  of  coupled  (s+d)-wave 
states.  A  brief  investigation  of  the  (s+d)-wave  state  on  a  square 
lattice  has  been  reported  previously  [26]  and  will  be  compared  with 
the  full  three-dimensional  results.  We  are  aware  that  the 
limitations  on  any  mean- field- theory  description  of  the  high-Tc 
materials,  namely  the  Brout  condition,  due  to  critical  fluctuations 
is  very  restrictive  [27];  however,  the  qualitative  agreement  of  the 
GL  thoery  with  experiment  deserves  mention. 

As  is  done  in  the  GL- theory  for  a  single  even-parity  order 
parameter,  we  write  the  free  energy  density  difference  between  the 
superconducting  state  and  the  normal  state  as  an  expansion  in  even 
powers  of  the  complex  gap  function  A(k) ,  which  is  related  to  the 
anomalous  thermal  average  <c£f  c  £.>  of  the  microscopic  theory  [28]A 
where  c^  is  the  electron  annihilation  operator  with  wave  vector  k 
and  spinet.  However,  for  the  multiple-order  parameter  case  we  must 
expand  A(k)  as  a  linear  combination  of  the  angular  momentum  basis 
functions  (Yj(ft)), 

2  2 
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consequence  of  this  is  to  induce  some  low- angular-momentum  s-d  2  2’ 

coupling  as  described  below.  Generating  the  invariant  terms  of  Xhe 
free-energy  density,  as  previously  described  [29],  we  can  write  the 
free-energy  difference  between  the  superconducting  and  normal  state 
for  a  tetragonal  lattice  as 
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Here  we  define  the  coherence  lengths,  £,,,  as  -  J<  /[2m.  ,|a.  j  ] , 

a  -  A, (T-T,),  where  j  refers  to  the  species  And  i  the  orientation 
(p  refers  to  the  xy-plane),  as  is  done  in  GL-theory  for  axial 
symmetry.  The  gauge- invariant  differential  operators  are^defined  as 
D  -  (3/3f  -  i^,  A  )  (f  -  x,y,z),  with  vector  potential  A  and  <f> . 

”  hc/(2e)  being  the  flux  quantum.  The  coupling  terms  in 
the  gradient  expressions  are  characterized  by  reciprocal  * ' 

H. ,  -  ty  /4m...  The  phase  angles  9,  and  9 


effective 
are  taken 
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the  pha^e  of  iJq,  thus  ensuring  the  gauge  invariance 
We  use  o/8n  to  represent  the  internal  magnetic  field 


masses,  .... 
relative  to ^9. 
of  Eq.  (2). 
energy  density 

Equation  (2)  has  been  subdivided  into  terms  arising  from  a  two- 
dimensional  analysis  of  the  square  xy-planes,  Jfs  and  J^g,  and  the 
additional  terms  required  to  analyze  tetragonll  systems,  and 
3^.  Reduction  of  the  symmetry  to  an  orthorhombic  point  group  adds 
an  additional  term  of  the  form  A_Aj,cos0^  to  Eq.  (2c)  and  destroys 
the  axial  symmetry  of  the  gradient  terms.  iIn  light  of  the  smallness 
of  the  orthorhombic  distortion  and  prevalence  of  twinning  in  the 
copper  oxide  superconductors  [30],  we  assume  the  AqA^  term  to  act  as 
a  perturbation  on  the  free  energy  of  the  tetragonal  lattice  and 
ignore  the  effect  of  reduced  symmetry  on  the  gradient  terms.  It  is 
interesting  to  note  that  the  two  d-wave  states  above  do  not  couple 
directly  to  each  other  up  to  order  i  -  2  of  the  relative  orbital 
angular  momentum  of  the  Cooper  pairs . 


Results 


We  have  performed  a  full  minimization  of  the  free  energy  with 
respect  to  the  A.'s,  9 s  and  the  vector  potential  A  to  obtain  a 
self-consistent  picture  of  the  thermodynamics  and  spatial  variation 
of  the  order  parameters  which  reproduces  the  dominant  features  of 
the  single-crystal  data  of  the  high-T  oxides.  Even  though  many 
parameters  appear  in  Eq.  (2),  we  understand  the  basic  physics  in 
simple  qualitative  terms.  The  simplest  scenario  is  that  of  the 


coexistence  of  a  highly  anisotropic  d  ^  2'state*  4^,  responsible 

for  the  quasi- two-dimensional  charac^er^of  these  materials,  with  a 
nearly  isotropic,  mixed  (s+d  2)-state,  possibly  characterizing  the 

"holon" -pair  hopping  wihtinzthe  RVB  picture  [31].  As  determined  by 
Kotliar  [11],  the  transition  temperature,  T^,  of  the  H-state  is 
higher  than  that  of  the  mixed  state.  A  schematic  picture  of  the 
relative  magnitudes  of  the  order  parameters  is  given  in  Fig.  1.  The 
relative  phases  are  0^  -  x/2  and  0j  “  n  near  the  transition 
temperatures.  The  small '“'amount  of  A  -  Aq  +  A2  state  persisting 
above  the  onset  temperature,  T,  is  a  consequence  of  the  small 
perturbation  to  Eq.  (2)  caused  by  a  shift  from  tetragonal  to 
orthorhombic  symmetry.  Perhaps  in  a  naive  way,  this  may  be  viewed 
as  adding  the  three-dimensional  character  necessary  for  the  onset  of 
superconductivity  [9],  The  existence  of  d-wave  states,  consequently 
gapless  superconductivity,  would  explain  the  large  observed 
Sommerfeld  coefficients,  while  the  multiple  transitions  of  these 
states  would  explain  the  two  specific  heat  anomalies  observed  near 
Tc  [32,33]. 

We  feel  that  the  peak  in  the  ultrasound  attenuation  results 
from  the  oscillatigns  of  the  relative  phases  0^  and  0 ^  about  their 
equilibrium  values  0^  -  ir/2  and  0«  -  it,  as  suggested  by  Kumar  and 
Wolfle  [13]  in  a  different  context.  Defining  w7  -  d  Jr^/ddi  (j  - 

1,2),  where  +  JL,  the  oscillation  frequencies  are  -*given 

by  L.  sq  1 

2  , 


i,  aa  ivuutaA  aiiv* 

:.  Defining  wf  -  d  Jr^/ae^  (j  - 
.llation  frequencies  are  J  given 


2  2  2 
wi  -  4A0A151 


u2  "  A0A2(A2  '  8A0A262  +  ^20A0  +  ^22A2^  (3b) 

There  will  be  a  sharp  onset  of  these  oscillations  at  T  which  will 

correspond  to  the  attenuation  peak  at  T  -  0.9  T  . 

We  next  consider  the  variation  of  the  upper  critical  field,  H  2 

with  orientation  and  temperature.  Using  a  straightforward 

variational  approach  on  the  linearized  form  of  Eq.  (2),  we  have 

derived  the  differential  GL  equations,  the  full  details  of  which 

will  be  presented  elsewhere.  For  the  sake  of  simplicity  we  assume  a 

(s+d  _)-wave  mixed  state  with  A„  -  A_  -  A  and  -  F  and 

i.  0  2m  Up  Zp  11m 

writ!  differential  equations  for  fields  11  parallel,  H*’ ,  and 

gerpendicular,  H  to  the  xy-plane.  For  H"  -  (H,0,0)  and 

A  -  ( 0 , - zH , 0 ] ,  we  have , 

A  -  (£  Hz)2  A  +  £2 (d2A  /dz2)  -  0  ,  (4a) 

m  ^mrinv  m  sm  m' 

A1  *  (*lp*invHz)2  A1  +  £izCd2V<iz2)  ”  0  '  (4b) 

1 

Similarly,  for  H  -  (0,0, H)  and  A  -  (0,xH,0),  we  have 

<o_- A„)A_  -  (a_£2+2Mno)(^__Hx)2A_  +  (a^2+2Mni )  (d2Am/dx2)  -  0, 


These  equations  are  decoupled  and  can  readily  be  solved  for  H 
within  the  harmonic  oscillator  approximation  to  yield  " 

4  -  •  4  -  <w!P>'1 

“c2  -  ■  Wf*i„v<4  +  2M02/V1'1  •  <5> 

Figure  2  gives  the  variation  of  the  critical  fields  with  temperature 
for  <  £  <  £.  .  For  the  upper  critical  field  is  always 

determined  by  ¥he  smaHest  coherence  length  £..  (0  K)  -  Jk.  For  H 
the  upper  critical  field  becomes  the  larger  between  H c2  and  H  2  as 
given  above.  This  may  explain  the  discrepancy  in  the  reported  cu  K 
values  of  the  in- plane  coherence  length  (?  (0)  -  34A,  €  (0)  - 

22A) ,  as  well  as  the  kink  in  the  H  „  data.  " 

The  variation  of  the  lower  critical  field,  Hc^>  with 
orientation  and  temperture  for  the2mixed  state  can  be  approximated 
by  the  expression  H  .  -  (®q/4w  A  £f)in  (*e£f)  [34],  which  is  valid 
for  large  values  of  tfie  GL-parametlr ,  ”  ^eff^eff’  For  c^^-s 

case  the  variation  of  the  internal  filla  occursmainly  in  a  region 
where  the  order  parameters  exhibit  their  maximum  values.  One  can 
therefore  obtain  the  penetration  depth,  by^casting  the  current 

relations  into  the  form  of  the  London  equation,  V  x  b  -  -  A  "  -A. 

The  results  for  H''^  and  H  ^  are,  6 


-2  -2 

2A  +  A,  +  A 
m  lz  2 


where  the  same  assumptions  on  Aq  and  A2  were  made  as  for  the 

calculation  of  H».  At  tempertaures  near  Tc  -  T^  the  lower  critical 

field  should  behave  as  A^  since  it  is  proportional  to  the  square  of 

the  order  parameter.  Consequently  the  anisotropy  of  should  go 

as  the  square  of  the  anisotropy  of  Hg*-  At  ^ower  temperatures  the 

influence  of  the  coupling  terms  A'  and  A  makes  predictions  more 

difficult.  The  anticipated  behavio?  of  H  ,  for  several  values  of 

the  coupling  terms  is  given  in  Fig.  3.  we  are  at  present  not  aware 

of  any  single-crystal  H  ,  studies  over  the  entire  temperature  range 

0  -  T  ci 

c ' 

Summary 

We  have  analyzed  the  thermodynamic,  magnetic  and  ultrasound 
attenuation  data  on  oriented  samples  of  the  high-Tc  superconductors 
within  the  context  of  anisotropic  Ginzburg- Landau  theory  for 
coupled,  even-parity  superconducting  states.  We  are  able  to  present 
a  consistent  interpretation  of  the  data  in  terms  of  the  coexistence 
of  a  quasi- two-dimensional  d-wave  state,  with  critical  temperature 
T-  -  T  and  a  more  isotropic  mixed  (s+d)-wave  state  with  critical 
tempertaure  T  <  T  .  We  predict  the  possibility  of  a  "kink"  in  the 
temperature  dependence  of  the  lower  critical  field  near  0*9Tc,  which 
should  be  tested  by  experiments  on  single  crystals. 
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Figure  1.  Schematic  temperature  dependence  of  the 
superconducting  order  parameters,  where  Am  is  for  the  mixed 
(s+d)-state  and  for  the  pure  d  ^  2'sCate<  anc*  T  -  a  - 

T  are  the  critical  temperatures  i>f  the  mixed  and  pure 
spates,  respectively,  and  T  -  b  -  T  is  the  onset  temperature. 

Figure  2 .  Schematic  temperature  dependence  of  the  upper 

critical  field,  H  _.  The  dashed  curves  are  not 

experimentally  observable.  H'L  is  the  field  parallel  to  the 

ab-plane,  and  _H  «  is  the  field  parallel  to  c-axis.  T  -  a  - 

T  and  T  -  b  -  T. 
m 

Figure  3 .  Schematic  temperature  dependence  of  the  lower 
critical  field,  H  ...  The  dashed-curves  represent  the  effect 
of  the  coupling  terms  A  and  A  T  -  a  -  T  and  T  -  b  -  T. 
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